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We study the scattering of magnetoinductive plane waves by internal (external) capacitive (in- 
ductive) defects coupled to a one-dimensional split-ring resonator array. We examine a number of 
simple defect configurations where Fano resonances occur, and study the behavior of the transmis- 
sion coefficient as a function of the controllable external parameters. We find that for embedded 
capacitive defects, the addition of a small amount of coupling to second neighbors, is necessary for 
the occurrence of Fano resonance. For external inductive defects, Fano resonances are commonplace, 
and they can be tuned by changing the relative orientation/distance between the defect and the SSR 
array. 
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I. INTRODUCTION 

Metamaterials (MMs) are novel artificial materials 
characterized for having negative dielectric permittiv- 
ity and negative magnetic permeability over a finite fre- 
quency range, endowing them with unusual electromag- 
netic wave propagation properties [TJ [2] ■ One of the most 
studied MMs is a metallic composite structure consisting 
of arrays of wires and split- ring resonators (SRRs). The 
theoretical treatment of such structures relies mainly on 
the effective-medium approximation where the compos- 
ite is treated as a homogeneous and isotropic medium, 
characterized by effective macroscopic parameters. The 
approach is valid, as long as the wavelength of the electro- 
magnetic field is much larger than the linear dimensions 
of the MM constituents. The magnetic properties of SRR 
arrays have been explored in a number of works [5Hl0j. 

Among the numerous aspects to explore with a new 
and promising electromagnetic material such as MMs, 
are its optical transport properties. For a SRR periodic 
array, it is only natural to explore the propagation of 
long-wavelength excitations, in the presence of one or 
few "impurities" that break its translational invariance. 
In the absence of any defects, the eigenmodes of a lin- 
ear periodic SRR array are termed "magnetoinductive 
waves". The presence of judiciously placed defects make 
possible interesting resonance phenomena, such as Fano 
resonances (FR), where there is total refiection of plane 
waves through the impurity region, in an otherwise pe- 
riodic potential. In a typical FR system, the wave prop- 
agation in the presence of a periodic scattering poten- 
tial is characterized by open and closed channels. The 
open channel guides the propagating waves as long as 
the eigenfrcquencies of closed channels do not match the 
spectrum of linear waves. The total reflection of waves 
in the open channel occurs when a localized state orig- 
inating from one of the closed channels resonates with 
the open channel spectrum|ll]. The FR effect can be 
used to control the optical response in novel, man-made 
materials, such as magnetic metamaterials. 

In this work, we consider Fano resonance effects due 



to a few defects whose position, with respect to the peri- 
odic SRR array can be easily tuned, making this type of 
configuration an interesting one to probe experimentally. 

This paper is organized as follows: Section II reviews 
the basic equations describing a one-dimensional, weakly- 
coupled SRR array in the absence of defects. In section 
III we examine a number of simple defect (s) configura- 
tions and obtain the transmission of plane waves across 
the impurity region, focusing on possible Fano resonance 
phenomena and, finally in section IV we summarize our 
findings. 



II. SPLIT-RING RESONATORS (SRR) ARRAY 

Let us consider a one-dimensional, periodic array of 
split-ring resonators (SRRs), in the absence of nonlin- 
earity, driving and dissipation. The most simple form of 
a split-ring resonator is that of a small, conducting ring 
with a slit. In general, each SRR unit in the array can 
be mapped to a resistor-inductor-capacitor (RLC) circuit 
featuring self-inductance L, ohmic resistance R, and ca- 
pacitance C built across the slit. In our case, we will 
consider the case of negligible resistance R. Thus, each 
SRR unit possesses a resonant frequency Wq ~ 1/ VLC. 
In the array, each SRR is coupled to their nearest neigh- 
bor via mutual inductionjS]. 

The evolution equation for the charge Qn residing on 
the nth ring is 



LQn + ^^nmQ,n + 7^ = (1) 



where Mnm is the mutual induction between rings n and 
m. We cast this equation in dimensionaless form by defin- 
ing T = iO(it, qn = Qn/CUo, Xnm = Mnm/L, whcrc Un is 
a characteristic voltage across the slits. 

The dimensionless evolution equation for the charge qn 
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FIG. 1: Two typical configurations for tlie SRR array. 



residing on the nth ring reads now 




q7i = 



(2) 



where Xnm denotes the ratio of the mutual inductances 
between the nth and the mth ring to the self inductance 
of the rings, and decreases as the inverse cube of the ring- 



to-ring distance, 



1-3 



(to 7^ n). In the limit 



of weak coupling (large distance between SRR units), it 
is customary to take A„„i = A Sn.m- 

The two most simple SRR configurations are shown in 
Fig.l. In one of them (top), all couplings are positive, 
while for the second (bottom), they are all negative, as 
a result of Lenz's law. The stationary modes of Eq.([2| 
are obtained by posing a solution in the form <7„(t) — 
g„exp(zf2t). This leads to the stationary equations. 



(3) 



Magnetoinductive plane waves of the form — 
A exp(ifcn), lead to the dispersion relation 
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1 + 2 J2m>0 '^Om COs(fcm) ' 



(4) 



We must impose the physical condition fi^ > 0. Using 
that Aom has the form Aom = X/m^, one obtains: 



1 2 
2C(3) < ^m)' 



(5) 



or, -0.41595 < A < 0.554605, where C(s) is the Riemann 
Zeta function ({s) = J2T=i '^^ ^^^^ s^^' hmits 

in condition ([5| change a bit when one assumes a weak 
coupling limit where few or even only one of the Aom is 
retained. 



III. MAGNETOINDUCTIVE DEFECTS AND 
FANO RESONANCES 

We consider now several simple cases involving defects 
coupled to the SRR. These defects are created for in- 
stance, by altering the electromagnetic characteristics of 
some rings, which can be achieved by changing one or 
more geometric features of the rings. Another way is by 




FIG. 2: Top: Capacitive defect embedded in SSR array, witli 
couplings to first-and second nearest neighbors. Bottom: (a), 

(b) and (c) show transmission coefficient vs wavector, for sev- 
eral parameter values: A = —0.1, A' — (1/8)A,A — 0.7 (a), 
A = 0.4, A' = 0.21, A = -0.5 (b), A = 0.4, A' = 0.21, A = 0.5 

(c) , A = 0.4, A' = 0.21, A = 1.0 (d). 



altering their relative coupling with respect to the array. 
In the case where nonlinearity and dissipation are not 
considered, two simple choices are to alter the slit width 
of a ring which alters its capacity only, affecting the value 
of its resonant frequency. Another simple choice is to 
couple the array to an external defect (s) consisting of a 
ring(s) placed outside the SRR array. This affects the 
value of the coupling between the SRR and the external 
ring(s) only. 



A. Single capacitance impurity 

The first case we study is a single capacitance impu- 
rity embedded in the array (Fig. 2), which without loss of 
generality we place at n = 0. This impurity is created 
by changing the slit width of the ring at n = 0. The 
stationary equation for this case is 



E 



^^{qn + Xnmqra) + (1 + (5„,oA)(7„ = 0, (6) 



where A is the change in the (dimensionless) resonance 
frequency at the defect position. When the slit width 
tends to zero, its capacitance diverges, making the res- 
onant frequency approach zero. This implies A — > — 1. 
On the other hand, when the slit width is large, the lo- 
cal capacitance goes to zero, and the resonant frequency 
diverges, implying that A — ^ oo. Thus, —1 < A < oo. 

Usually, embedded defects coupled locally do not lead 
to Fano Resonance phenomena (perfect plane-wave re- 
fiection). In the SRR array however, the couplings are 
dipolar and therefore, long-range, and non-local effects 
have to be considered. For computational simplicity 



we will work with couplings to first-and second nearest 
neighbors only. The hope is that interlerence between 
the path through nearest neighbors and the one through 
next-nearest neighbors will give rise to Fano Resonance. 
Equation Q becomes: 



[qn + HQu+I + Qn-l) 



(l + A„o)gn = 0. (7) 



where A' is the coupling among next-nearest neighbors. 
We pose a plane wave solution ol the form 



qn 



Re 



-ikn 



n < -1 
n > 1 



(8) 



After replacing this ansatz into Q, one obtains after a 
little algebra: 
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l-|-2Acos(fc)-h2A'cos(2A;) 



and 



where, 



t EE \T/I\' 



A + AB 



A + AC 



(9) 



(10) 



A = -2isin(A:)[2A'^cos(3/c) + 6AA'(Acos(fc) + A'cos(2fc)) 

+A(A2 + 3A'2)] 
B = -2isin(A:)[ 2AA'(Acos(fc) + A'cos(2fc)) + AA' ] 
C = [ l + 2Acos(/c) + 2A'cos(2fc))(A2 + 2AA'e-*'= 

-A'2 + 2A'^ cos(2fc) ] 



(11) 



Now, for the SRR case where dipolar interactions fall 
with the inverse cube of the distance between SRR units, 
we have A' = (1/8) A. After inserting this into Eqs.([9]) 
and ( 10 ) and after imposing f]^ > 0, we conclude that 



-(4/9) < A < (4/7) is the relevant coupling regime for 
this SRR configuration. Inside this regime it is possible 
to have a relatively weak Fano resonance, as shown in 
Fig. 2a. As expected, when coupling to second nearest 
neighbors is neglected (A' = 0), there is no FR at all, 
since in that case the transmission becomes 



t = 



(2Asin(/c))2 



(2Asin(fc))2 + A2(l + 2Acos(fc))2 



(12) 



which can only be zero at k — 0, n. For more general 
A, A' values, it is possible to have zero, one, and even two 
different strong Fano resonances, as shown in Fig. 2b, 2c 
and 2d, respectively. 



B. Single inductive impurity 

Another simple case is a single inductive impurity em- 
bedded in the array. This is achieved by placing a SRR 
unit at different distances from the left and right por- 
tions of the array, leading to asymmetrical couplings Xl 




FIG. 3: (Color online) Top; SRR array with external induc- 
tive defect. Bottom: Transmission coefficient vs wavevector, 
for different values of the relative coupling defect: Ai/A = 
1.25(circles), 1.5(squares), 1.75 (rhombi) and 2.0(stars). 



and Afl. Only couplings between nearest-neighbors are 
asumed. This case was already considered in ref|12). In 
this case, it can be proven (not shown) that there is no 
FR for any choice of A^, A^; the transmission vs wavevec- 
tor curve shows a single maximum at fc = 7r/2. For the 
special symmetric case Xl — A^, there is a transmission 
resonance {t — 1) at k = Tr/2. 

We move now to simple cases where the defect (s) lie(s) 
outside the SRR array, and thus, their coupling can be 
tuned by changing their distance to the array. 



C. Single external inductive defect: First case 

The first case of this type is that of a single external in- 
ductive defect coupled to the SRR at n = via coupling 
to first and second nearest neighbors via coupling param- 
eters Ai and A2, respectively (Fig. 3). From the geometry 
of the configuration, it is easy to obtain the relationship 



A2 
A 



2/3 



-3/2 



(13) 



For simplicity, we assume only nearest neighbor cou- 
pling among the SRR units. For this approximation 
to be consistent, one must impose Ai > (1/8) A and 
also A2 > (1/8)A. From Eq.([l3]), this implies Ai/A > 
Max{(l/8), 3-3/2} ^ 3-3/2 ^ 0.192. 

The stationary equations for this case are 



- ^'^[Hqn+i + qn-i) + qn 

+qe{XlSn,0 + X2Sn,±l)] + qn = 

-fl'^{Xiqo + X2q±i+qe) + qe = 



(14) 
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where Qe is the charge residing on the external defect 
ring. We assume a plane wave solution of the form 



Re 



-ikn ^ < Q 

n> 1 



(15) 



The transmission coefficient t 
uration is 



V A , 



4 ((»)(ij 



1+ ¥ 



|T//|2 for this config- 



cos(fc)^ sin(fc)) 



lf) + (l 



2isin(2fc) 
(16) 

where, A2 is given by Eq.(13) for the SRR system. The 



most interesting feature of this transmission coefficient 
is the presence of a Fano resonance, whose position 
varies with the value of Ai/A. The FR occurs when 
A1A2 + (A^ + A2)cos(fc) — 0, which is possible from 
Eq.([l3|, for all < |Ai/A| < 2.449. Therefore, the rel- 
evant interval where FR exists in our model is given by 
0.192 < Ai/A < 2.449. This can be achieved by simply 
changing the distance between the SRR array and the 
external defect. Note from Eq.([l6| that, if we neglect 
the effect of A2 , the position of the FR remains locked at 
fc = 7r/2. 

To get a better feeling as to what would be expected 
to see in an actual experiment for this system, we per- 
form a numerical simulation of a wide pulse impinging on 
the inductive impurity region. We used an array of 10'^ 
rings and resort to a symplectic algorithm to trace the 
time evolution of an initial broad gaussian pulse m„ (0) = 
Aexp{—a{n — no)^) exp(ifc(n — no)), where uq = —300, 
A — 0.1 and a — 0.001 which means a width of about 
120 sites, in order to simulate a plane wave with a well- 
defined k. The rings coupled to the inductive defect are 
placed at n = 0,±1. Figure [3] shows numerical simula- 
tion results for the transmission coefficient vs wavevector 
for several A'/A values, with A = 0.4. The agreement be- 
tween analytics (Eq.(16)) and numerics is excellent, and 



the Fano resonances are clearly shown. The numerical 
discrepancies at wavevectors close to A: = 0, tt are due 
to the long integration times needed since the pulse is 
quite slow at these k values. In Fig. 4 we show the out- 
put profile of the gaussian pulse for different fc-values, 
after an integration time t — 1500. The transmission in- 
creases monotonically up to fc ~ 0.389 7r(0.868 tt) from 
the left (right), then decreases steadily until it vanishes 
completely at k « O.GOItt. 



D. Single external inductive defect: Second case 

The next case we consider is a variation of the previous 
case, where the external defect is now located halfway be- 
tween two SRR units (Fig. 5). Without loss of generality, 
we take the external defect coupled symmetrically to the 
units at n = and n — 1, with coupling A'. As before, 
the units in the array interact through nearest-neighbor 
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FIG. 4: (Color online) SRR array with external inductive 
defect: Output profile after t — 1500 as a function of wavector 
and ring site, for Ai/A — 1. Impurity site is located at n = 0. 



couplings A only. In order for this approximation to be 
consistent, and taken into account the dipolar nature of 
the inductive couplings, we need |A'| > (1/8)A. 
The stationary equations read 

- n'^{X{qn+l + qn-l) +qn+ A'(Je('5„,0 + <5„,i)) + Qn = 

^n\\'{qo + qi))+qe = (17) 

where Qe is the charge on the external defect ring. We 
pose a plane wave solution of the type 



R e 



— ikr, 



n<~l 
n>2 



leading to a transmission coefficient: 

(e''=-l)(A'2 + 2A2 cos(fc)) 



t ; 



T 


2 
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2(A'2 + A2 cos(fc)(l 



)) 



(18) 



(19) 



Fano resonances are possible when A'^ -I- 2X^ cos(fc) = 0. 
This implies, kp — — (1/2)(A'/A)2, which is possible only 
if |A'/A| < a/2- This constraint plus the consistency con- 
dition, give us the relevant coupling parameter window 
for this SRR configuration: 



(1/8) < 



< 



(20) 



Figure 5 shows some transmission curves for parameter 
values inside and outside this window. The position and 
width of the Fano resonance depend on the ratio |A'/A|, 
which can be externally tuned by changing the distance 
between the defect and the array. 

Figure 6 shows output intensity profile after propaga- 
tion oi t = 1500 as a function of wavevector and site 
position. 
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FIG. 5: (Color online) Top: External defect coupled sym- 
metrically to the SSR array. Only couplings to first nearest 
neighbors are considered. Bottom: Transmission vs wavevec- 
tor for several X' /X values: 0.5 (circles), 1.0 (rhombi) and 1.5 
(triangles) . 
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FIG. 6: (Color online) Fano resonance: External defect cou- 
pled symmetrically to two array units: Output profile af- 
ter t = 1500 as a function of wavevector and ring site, for 
Ai/A= 1. 



E. Two external inductive defects 

The final array we consider consists of two external 
defects, each coupled to a single unit of a SRR array. 
The units of the array are coupled to first neighbors only 
via coupling A. Between both defects, there are L units 
(Figj?]). Without loss of generality we assume the defects 
to be coupled to units at n = and n = L, with couplings 
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FIG. 7: Two external defects separated by L units and cou- 
pled to the SSR array. Transmission vs wavevector for several 
separations L: L = 2(a), L = 3(b), L — 5(c) and L = 20 (d) 
(in all cases A = Ai = A2 = 0.4). 



Ao and A^, respectively. The stationary equations read 

- ri^[A((7„+i + qn-l) + Qn + AoQo^n.O 



+A 



-f72(Aogo + Qo)+Qo = 

-n''iXLqL + QL) + QL = (21) 

where Qo and Q l are the charges residing at the external 
rings coupled to units n = and n = L oi the array. As 
before, we pose a plane wave solution 



/ e''^" + R e-**^" n < 
T e**^" n> L 



(22) 



The transmission coefficient is given in this case by 



[\l + 2zA- 



A^(e^'fe 
'sin(2A;)][A' 



1) 

2iA2 sin(2/c] 



_X2^2p2«fe(L+2) 



(23) 

which is symmetric under the exchange between the ex- 
ternal rings, depending only on |Ao/A|, |Al/A| and L. 
We see from Eq.(23) that, regardless of the separation 



L between the external defects, there is only a single 
Fano resonance located at k ~ n/2. Figure 7 shows ex- 
amples of transmission vs wavevector plots for different 
defects separation. However, it should be noted that, 
in a more realistic setting that considers the presence 
of dissipation effects, one would expect to observe some 
broadening of the resonances, as well as a decrease in 
their heights, as the separation between defects increases. 
For SRRs coupled to nearest-neighbors only, this effect 
is more pronounced at the band edges, and minimized at 
mid-band I2J. For couplings beyond nearest-neighbors, 
the situation is more complex 13 . 
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IV. SUMMARY 

In conclusion, we have examined Fano resonance ef- 
fects in a SRR array coupled to internal (capacitive) and 
external (inductive) defects. In the case of internal or 
embedded defects, the presence of coupling to second 
neighbors, in addition to the customary first neighbor 
coupling, is necessary to induce the FR phenomenon. For 
external defects, the FR phenomenon depends mainly of 



the geometric configuration between the defect (s) and 
the SRR array. The position and strength of these FR 
could be easily tuned in current experiments, and con- 
stitute a clear example of the possibility of controlling 
the transport of electromagnetic waves across a magnetic 
metamaterial. 
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